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The introduction of intermediate states in the dynamics of the voter model modifies the ordering 
process and restores an effective surface tension. The logarithmic coarsening of the conventional 
", voter model in two dimensions is eliminated in favour of an algebraic decay of the density of in- 

f- ») ■ terfaces with time, compatible with Model A dynamics at low temperatures. This phenomenon 

' is addressed by deriving Langevin equations for the dynamics of appropriately defined continuous 

£\j . fields. These equations are analyzed using field theoretical arguments and by means of a recently 

proposed numerical technique for the integration of stochastic equations with multiplicative noise. 
We find good agreement with lattice simulations of the microscopic model. 
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I. INTRODUCTION 



o ; 

The theory of phase ordering dynamics is a prominent tool to study and classify the non-equilibrium behaviour of 
interacting particle systems Applications can be found in very different fields, ranging from the study of binary 
! 1 1 , mixtures [2] over reaction-diffusion systems Q to models of social dynamics and opinion spreading Q • 

In many physical systems, phase ordering occurs through defect annihilation and domain coarsening, as it can 
be observed for example in the Ising model after a quench from the disordered phase into the ordered one at low 
,i_S i temperatures Here, the typical length scale L(t) of clusters of aligned spins grows algebraically with time, following 
the general behaviour L{t) ~ i 1 / 2 , where e.g. z = 2 for non-conserved dynamics. The numerical value of the dynamical 
exponent z for any given model here depends on universal features of the model, such as symmetries and conservation 
laws, and on the relevant external parameters. 

The low-temperature ordering in Ising-like models exhibits another very general property: the existence of surface 
tension between domains, favouring the formation of smooth interfaces minimizing the curvature. This motivates the 
I notion of curvature-driven coarsening. In presence of smooth interfaces, we also expect the density of interfaces n(t) 
to be inversely proportional to the domain length n(t) ~ L^ 1 , thus decaying with similar algebraic law n(t) ~ t~ s 
I with S = l/z. 

Other systems exist in which coarsening takes place without surface tension [5[; the growth of clusters is there 
£-««. ' purely diffusive, driven by interfacial noise. A typical example is provided by the voter model ((J, in which spins can 
7-H , take one of two possible states, and Hip with a probability linearly proportional to the fraction of unlike neighbors. 
In dimension d = 1, the voter dynamics is the same as zero-temperature Glauber dynamics of the Ising model [7], 
but increasing space dimensionality the ordering slows down and becomes logarithmic in d = 2 (n(t) ~ l/log(f)), 
\Q ■ whereas an infinite system does not order for d > d c = 2. The notion of logarithmic coarsening in two dimensions 
here refers to the behaviour of the density of interfaces with time, whereas the dynamical exponent, relating to the 
correlation length and properties of the dynamical structure factor, is still found to be z = 2 This is due to the 
critical nature of the voter model, as unveiled using field theoretic approaches 0, [Tfj| . The voter dynamics corresponds 
^ , to a massless field theory without field renormalization (from which z = 2), but it admits renormalization of other 
quantities, including the density of interfaces n(t). Exact renormalization group calculations predict a behavior that 
agrees perfectly with both analytical solutions and numerical simulations [111 . 1X21 ] . 

Recent research in opinion dynamics and sociophysics has amplified the interest in these two distinct classes of 
non-equilibrium systems, and has prompted their use in the study of the collective behavior of social individuals (4j. 
For instance, voter- like models with multiple states have been used to model groups of political parties [H, [HI, the 
spreading of linguistic conventions [TBI . [l6l. llTj , while the introduction of counters fl8j or inertia in the spin dynamics 
[19( have been used to describe the effect of learning and memory in social interactions. Interestingly, the mere 
introduction of an intermediate state, without any major modification of the voter dynamics itself, has been seen to 
be sufficient to avoid logarithmic coarsening in two dimensions, restoring an algebraic decay of the density of interfaces 
and an effective surface tension [l5|, [l8[. As we will show below, the qualitative asymptotic behavior does not change 
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TABLE I: Transition rates of the two models considered. The 2-state (spin 1/2) model is the usual voter model, while the 
3-state model (spin 1) is the one with an intermediate state. The quantities f±,fo are the fractions of ±,0 spins among the 4 
nearest neighbours of a given spin. 

if one allows for more than one intermediate state, provided that only two absorbing states are present and that they 
correspond to the extreme states. 

A better understanding of this intriguing phenomenological picture can be obtained by regarding the voter model as 
a special point of a larger class of non-equilibrium systems with two absorbing states, exhibiting a Z2-symmetry, see 
SUS I1H- Starting from this viewpoint, Al Hammal et al. [22j have derived a simple phenomenological field theory, 
expressed in terms of a non-linear Langevin equation, which captures the essential characteristics of the dynamical 
behavior of these models, and which faithfully reproduces the logarithmic coarsening dynamics of the voter model. 

The aim of our work is to provide a theoretical explanation for the restoration of curvature-driven coarsening by 
intermediate states, and to corroborate the microscopic view with a continuum approach able to clarify the link with 
known results on models in the the so-called generalized voter (GV) class 0, (22[. In the following we show analytically 
that voter models with intermediate states (this includes e.g. Naming Games [TEl [l6j) can be regarded as members 
of this class, indirectly providing an example of non-trivial microscopic models belonging to the GV class. We will 
derive phenomenological equations for the 3-state voter model, and describe how such equations are, within certain 
limits, useful to describe the dynamics of models with a larger number of intermediate states as well. These Langevin 
equations are studied both numerically and more formally using field theoretical arguments. We show that despite 
the presence of interfacial noise, in two dimensions the asymptotic behavior is compatible with the low (but non- 
zero) temperature regime of Model A relaxational dynamics. We also show that, unlike in the Ising model, interface 
roughening is not expected to occur in higher dimensions. 

II. VOTER MODELS WITH AND WITHOUT INTERMEDIATE STATES 

In the voter model Q spins are arranged on a network with one spin at each node, and each spin can assume two 
possible values, e.g. ±1. At each step of the microscopic dynamics one spin is selected at random, and subsequently 
assumes the state of one randomly chosen nearest neighbour. We introduce an intermediate 0-state, maintaining a 
dynamics similar to that of the voter model. One spin s and one of its neighbours s' are chosen randomly at each 
time step. If the spin selected is in state s = +1 and the neighbouring one in state s' — — 1, then s assumes the value 
s = 0. Similarly for s = — 1, s' = 1 the first spin will be set to s = 0. If s = 1 and s' = 0, then s will assume the state 
s = with probability 1/2 and remain at s = 1 otherwise. Similarly for s = —l,s' = 0. If s = 0, then s will adopt 
the state s' if s' = ±1, and will flip to either s = 1 or s = — 1 with equal probability if s' = 0. The exact dynamical 
iteration prescriptions and transition rates for the voter model and the modified 3-state model are summarised in 
Table HI Note that with these rules intermediate states do not persist in time, and the only possible absorbing states 
are those with either all spins up or all spins down. In this work we will focus on models defined on d-dimensional 
regular lattices of lateral size L, and most of the numerical analysis is carried out in the two-dimensional case, where 
we impose periodic boundary conditions. Interaction occurs with the set of the four nearest neighbours of any given 
spin. 

The 2-state and 3-state models behave in the same way in one dimension, since the presence of intermediate 
states only has the effect of renormalizing the width of the domain walls by a finite factor (see [H, [l8| for similar 
considerations). In two dimensions the behaviour of the model with intermediate states is instead very different from 
that of the standard voter model. This can be seen in Fig. [1] where we report the temporal evolution of the density 
of intermediate states for the 3-state model. We note here that spins in the intermediate state are found at the 
interfaces of domains of up and down spins respectively, so that the density of intermediate states is a measure of the 
density of interfaces (or domain boundaries). In the standard two-dimensional voter model, the interface density is 
known to decrease logarithmically, while the modified (3-state) model shows algebraic decay, i.e. n(t) ~ t~ s , where 
6 has been reported in the range of 0.45 to 0.5 [l5j]. The correlation length £, which we have computed from the 
two-point correlation function, grows with a similar algebraic law (see inset of Fig. [I}. Both measures seem to suggest 
a dynamical exponent z = 2, even if from these numerical values the contribution of logarithmic corrections cannot 
be excluded. These were also the conclusion of previous studies [HI, EH) an d are corroborated by the marginality of 
voter dynamics in two dimensions. The point will be clarified in the next sections by means of a more theoretical 
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FIG. 1: (Colour on-line) Density of intermediate states (main panel) for the 3-state model. Results from one sample of a direct 
microscopic simulation of a system of size 2000 2 are shown, along with measurements from numerical integration of Eqs. 
where we have ignored the noise terms (see text; system size is 500 2 , averages over 10 runs with independent random initial 
conditions), and from Eq. |(9]| at a = 1/2, D — 1 and different noise strengths 7 (curves at 7 > are obtained from systems 
of size 500 2 , integration method as explained in the main text, averages over 10 samples are taken; results for 7 = are from 
Euler-forward integration, time-stepping of 0.05, 5 samples of size 1000 2 ). The inset shows the behaviour of the correlation 
length £ of the field (f> (obtained from the circularly averaged correlation function). All curves have been rescaled by arbitrary 
constant factors for optical convenience. One unit of time corresponds to one update per spin on average. 
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FIG. 2: Snapshots of the time evolution of the model with intermediate states. Sites i with Si — —1, 0, 1 correspond to black, 
white and grey respectively. Data is from simulations of the microscopic model, only a region of N = 100 2 sites is shown. The 
run in the upper row started from random initial conditions Sj = ±1, snapshots are taken after 100,200,400 and 800 sweeps 
over the system. Lower row shows evolution of a droplet (snapshots are taken at time t = and then after 500, 1000 and 2000 
sweeps) . 



analysis. 

Fig. [5] displays snapshots of the coarsening dynamics started from a disordered configuration (top), and during a 
'droplet experiment' in which a circular domain of up spins is left to evolve in a sea of down spins (bottom). 
In the pictures it is visible that despite the voter-like update rule, the 3-state model develops an effective surface 
tension, as already noted for similar dynamics with intermediate states, such as the Naming Game model [l6j . the 
noise-reduced voter model [l8[ , and models of bilingualism HU . 

To investigate the role of the interfacial noise clearly visible in the snapshots in Fig. [2 we have analyzed the roughness 
properties of the interface between two domains. Let us consider a flat interface between +1 and —1 spins and study 
its evolution in time under the dynamics of the model (see 2J| for details of such experiments in the context of other 
models). Fig. [3] displays the behaviour of the average fluctuations of the width of the resulting interface . More 
precisely, if we call hi the displacement in the direction perpendicular to the interface at position i (the latter is 

1/2 



measured along the interface), the width is defined as W — 



where L is the linear 
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FIG. 3: (Colour on-line) Data collapse for the interface width w(L, t). The time is rescaled by L z (having assumed z = 2) and 
the width is rescaled by L a with a — 0.5. This generates an initial power law growth t 13 with /3 ps 0.25 (dashed line). This 
rescaling is thus consistent with the Edwards- Wilkinson universality class of surface growth. 



size of the interface. As the dynamics progresses W assumes a stationary value asymptotically. Considering systems 
of different sizes, we find that this saturation value W does not remain finite in the infinite systems, but that instead 
the plateau value depends on the linear size of the system L. This indicates the presence of interface roughening 
phenomena. As in usual surface growth [24], the time evolution of the surface width exhibits two stages separated 
by a 'crossover' time t x . The initial growth follows a power law, W(L,t) ~ i' 3 , with growth exponent /?; then for 
t t x the width saturates at W sa t(L) ~ L a where a is the roughness exponent. Note that in the scaling theory of 
surface growth these exponents are related to the dynamic one by z = a/ (3. Figure [3] verifies the Family- Vicsek scaling 
relation W(L, t) ~ L a f (t/L z ), with a — 0.5 and z = 2. Accordingly we find f(u) ~ u@ with (3 w 0.25 at small times. 
The good data collapse generated by rescaling, suggests that the interface dynamics in our model is compatible with 
the Edwards- Wilkinson (EW) universality class of surface growth [24| . 

In summary, the numerical results indicate that the introduction of an intermediate state makes the voter dynamics 
similar to the low-temperature (T > 0) scaling of model A [25| . both as far as the domain growth and interface 
roughening arc concerned. 

Further generalisation to multiple intermediate states has already been considered. For instance, in the Naming 
Game [HI a large number of different intermediate states can be created. These are however all equivalent with no 
ordering or hierarchy among them, therefore the effective behavior of the system is expected to be the same as for 
3-state model. More interestingly we can define a voter model with non-equivalent intermediate states: let us consider 
S > 2 states in total, two of them being extremal absorbing states, and the remaining S — 2 ones being intermediate 
states. To implement this generalisation let us label the states by k = 1, . . . , S, with k = 1 and k = S being the 
extremal ones. As before the model is considered on a two-dimensional lattice, and at each time-step one spin and 
one of its four nearest neighbours are chosen at random. Let s G {1, . . . ,S} be the state of the chosen spin, and s' 
that of the neighbour. Then set p(s') — (s' — l)/(S — 1). As a consequence of the interaction, spin in state s is then 
set into state s + 1 with probability p(s')Q(S — s), and into state s — 1 with probability (1 — p(s'))Q(s — 1). Q(x) is 
here the step function with Q(x) = 1 for x > and Q(x) — for x < 0. It ensures that spins are never modified to 
take values outside the set {1, . . . , S}. As a consequence of the definition of p(s'), linearly increasing from to 1 as 
s' increases from I to S, the extremal states 1 and S have maximally polarising effects, whereas intermediate states 
induce flips of neighbouring spins only at a probabilistic rate. For S — 2 the model reduces to the voter model, while 
for S — 3 we recover the 3-state model as discussed above. 

Results from simulations of this generalised model for several values of S are reported in Fig. [4l while some 
snapshots of the dynamics for the model at S = 5 are shown in Fig. [5] As seen in the figures, the coarsening of 
the model is curvature-driven also when multiple intermediate states are present, and the density of intermediate 
states decays with an exponent close to the one already found for the 3-state model almost independently of the 
number of allowed states S. In simulations with many intermediate states (S of the order of ten or more), interfaces 
can be expected to have a finite size (typically O(S)). The definition of the measured quantity (density of spins in 
intermediate states) can hence be affected as the number of intermediate states is systematically increased, potentially 
limiting the accuracy of our numerical results and preventing one from a reliable measurement of the scaling exponent 
in simulations, even at comparably large lattice sizes (e.g. 1000 2 ). Still, results reported in Fig. 2] are consistent with 
algebraic coarsening at an exponent of 0.45 to 0.5. 

Other definitions of multi-state models are of course possible. For instance one could consider a model with S 
states (again with an internal ordering between them) with two extremal absorbing states (1 and S), but with 
uniform transition probabilities among neighboring states (e.g. a given spin in state s may change to state s + 1 
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FIG. 4: (Colour on-line) Coarsening of the model with multiple intermediate states. The plot shows the density of spins in 
intermediate states for the model at S — 3, 4, 14 (solid lines from bottom to top). The dashed line is n ~ t -0 ' 5 and is shown 
as an optical guide only. Each curve represents data from a single run of the microscopic dynamics on a lattice of 1000 2 sites. 
One unit of time corresponds to one update per spin on average. 

upon interaction with a spin in any state s' > s, with a rate which is independent of the specific state s'). We 
have tested such models (results not shown here) and have found algebraic domain growth with a coarsening 
exponent of about 0.45, so that we conclude that the precise definition of the transition rates between internal 
states is immaterial as far as the growth of the characteristic length scale of the resulting domain structure is concerned. 



III. CONTINUUM DESCRIPTION IN THE SCALING REGIME 

In order to complement simulation results and to obtain analytical insight into the dynamics of the voter model 
with intermediate states we develop a description of the 3-state model in terms of continuous field variables. As the 
microscopic model is neither diluted nor endowed with mobility, a standard derivation of continuous equations for 
coarse-grained degrees of freedom (2|| [53] is not straightforward. For the same reasons and because of the fact that 
the model considered here is effectively a spin-1 model, an exact microscopic mapping to field theoretic methods 
using second quantization, e.g. following the lines of [H, [2^], is difficult as well. 

We will therefore proceed to derive approximate stochastic evolution equations directly from the microscopic rates, 
based on the assumption that fields are Gaussian random variables. Because of the nature of the derivation, these 
stochastic equations do not describe coarse-grained degrees of freedom, but microscopic ones and differ from those 
obtained by adding diffusion to mean-field like terms [26| . The macroscopic description becomes clear when the slow 
modes are identified, by applying power counting and neglecting higher-order terms in the momentum space. This 
approach is somewhat similar to the derivation of effective field theories in condensed matter spin models. 

We start by considering a version of the model, defined as before on a 2-dimensional regular lattice, but now 
endowed with Q spins per site, and in which each spin can take any of the three states {—1, 0, 1}. The dynamics at 
each time step then proceeds by selecting one site and one of the £1 spins at that site at random, and then executing 
the interaction with a randomly chosen spin at a randomly selected nearest neighbour site. The £1 copies of the system 
hence become cross-correlated. 

In the limit f2 — > oo deterministic equations for the evolution of average local quantities such as the average local 
magnetization ((sj)) and the average local density of 0-states (1 — (sf )), can be easily obtained with the same rate 
equation approach applied by Krapivsky in the case of pure voter model [H], Averages (...) are here understood 
as averages over the copies of the microscopic dynamics, i.e. (si) is the average magnetisation at site i resulting 
from the Q spins at that site. In the limit of an infinite number of particles per site we replace (si) — > <f>(x,t) 
and 1 — (sf) — > tfj(x,t). Using such continuous variables is appealing but dangerous as it may introduce ultraviolet 
divergencies as well as a certain level of arbitrariness in the definition of the continuum limit of local operators (like the 
lattice laplacian) . During the course of the calculations presented in this paper we will not consider a continuum limit 
in space, and the variable x is assumed to have the granularity of an arbitrary lattice spacing as in standard lattice 
field theories. We expect that the spatial continuum limit could be in principle considered appropriately defining a 
diffusion constant D. 

If only a finite number Q of spins per site is considered, then progress can be made starting from the master 
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equation describing this ensemble of interacting copies of the system. Let us label the state of spin i in copy k 
of the system by s\ (t) at time t. We will denote averages over the Q. copies by cj)(x, t) in the following, i.e. we 

,!) (k). 



will use the substitution f2 _1 J2k=i s i '(tf ~> 4 > { x ->'t)- Note that while overbars indicate averages in a system with 
an infinite number of spins per site, averages denoted by a tilde are over a finite number f2 of copies only. It is 
then easy to verify that a change in the state of spin i in one of the Q samples corresponds to a transition of the 
type, <f>(x,t) — > (f)(x,t) ± and %p(x,t) — > ip(x,t) ± ^ in the quantities describing the average over the samples. 

This leads to a master equation describing the evolution of the functional probability distribution V[{(/)(x, t), ij)(x, *)}]• 
Expanding to next-to-leading order in fi -1 , we get the following multivariate functional Fokker-Planck equation for 
the functional probability distribution V 
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where we have re-scaled time by a factor of il 1 , i.e. t — > t/fi. V 2 is here the lattice Laplacian, i.e. W 2 (f>(x,t) = 
3 — < / > ( :c j' : )}: where the sum over y extends over the four nearest neighbours of x. Noise terms in the 

corresponding Langevin equations for <j) and ip are hence of order fi -1 / 2 . They reflect the fact that averages in a 
model with only a finite number f2 of spins per site still remain stochastic quantities. Stochasticity only vanishes if 
£1 — > oo, and deterministic mean- field equations for <j> and ip are obtained (these would be Eqs. @ below, with the 
noise terms removed). Our approximation consists in describing the model with only one spin per site (i.e. the model 
introduced and discussed in the earlier sections) by the following two coupled Langevin equations 



f t 4>(x, t) = \<j>{x, t)i>{x, t) + i±« V 2 4>(x, t) + n(x, t), 

^(x, t) = -ty(x, t) + ^1 V 2 ^^, t) + £(x, t), 

where r)(x,t) and £(x,t) are Gaussian noise terms, with correlation given by 

(t){x, t)ri(x',f)) = (^(x, m(x',t')) = \ {i>(x, t)+l- 4>{x, tf - fa, t)V 2 0(x, t)} 5(x - x')6(t - t'), 

(n(x, t)Z(x',t')) = (C(x, t)r)(x' , t')) = | {-3ip(x, t)4>{x, t) + (1 - 3^(x, t)) W 2 (f>(x, t) } S(x - x')5(t - t'). 



(2) 



(3) 



We have therefore effectively disregarded higher-order terms in the expansion in f2 -1 , even if the model actually 
corresponds to = 1. Nevertheless, as we will show below, Eqs. ^\ -[3]) turn out to be a faithful mapping of the 
original dynamics of the 3-state model onto a lattice theory with continuous fields at each grid point, capturing the 
essential feature of the original model. It is here worth stressing again that -similar to the voter model- our microscopic 
dynamics do not allow for microscopic mobility of agents (e.g. hopping or interchange of neighbouring particles), thus 
the diffusion-like terms in @ are non-trivial features emergent from the microscopic model. In particular the diffusive 
terms in the deterministic parts of © are quadratic in the fields (e.g. ip"V 2 (f> and 0V 2 </>). Similar structures have been 
found in the context of other spatial systems for example in (30j . 

A study of a version of our model including microscopic mobility can be performed along the lines of (3l| , here local 
mixing is ensured by scaling the hopping rate suitably with the system size, and an effective description in terms of 
continuous fields can be obtained by means of an expansion in the inverse system size. This leads to equations which 
are similar to but not identical with Eqs. ([2]) and ([3]) ,32|. More specifically one finds noise amplitudes of the form 
© with the terms proportional to V 2 removed and with an additional overall prefactor TV -1 (N = L 2 being the 
system size) . The deterministic terms take a form very similar to those reported in @ , but with decoupled diffusive 
terms (oc V 2 (f>, V 2 ip) instead of objects of the form ipV 2 (f> and <\N 2 <\>. 

If we were to repeat the same procedure as above for the case of the pure voter model with two possible states (and 
without mobility or particle interchange), the resulting equation would be 



d_ 



(f>(x,t) = \J 2 (j)(x,t) +T](x,t), 



(4) 



where (rj(x,t)r](x' ,<')) = {l — <t>{x,t) 2 — 2(f>(x,t)V 2 (f)(x,t)} 5(x — x')S(t — t'). Apart from the diffusion term in the 
noise, Eq. (H|) is identical to the field-theoretic equation proposed by Dickman et al. for the voter model [33| . Naive 
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power counting on dimensional grounds reveals that the diffusive terms involving Laplacians in the noise correlator 
are irrelevant in any dimension in the long wavelength limit. Long-time large-scale properties of a microscopic voter 
model are therefore well described by the effective Langevin equation 

-cj>(x, t) = Z?V 2 <Kx, t) + V 1 - <f>(x, t)%(x, t), (5) 

with (rj(x,t)ri(x' , i')) = 5(x ~ x')5(t — t'). Note that the same equation could have been obtained also by adding 
the diffusive term to a zero-dimensional mean- field Langevin structure [331 ] - Despite its validity, such coarse-grained 
derivation misses a clear microscopic justification. 

Similarly, on the base of dimensional analysis we could neglect the terms proportional to V 2 ^ in Eq. ([3]), providing 
slightly simplified noise terms with the same absorbing states. 

As we will see, numerical integration of Eqs ([2]) and ([3]) is challenging (simple forward integration in discrete time 
steps renders the field-dependent variances of the noise terms negative with finite probability), therefore it is useful 
to further reduce their complexity by means of some physically meaningful approximation. 

Provided we focus on the asymptotic dynamics of the system, it is possible to reduce Eqs. and to a single 
equation. We here exploit the observation that intermediate states cannot proliferate. Due to their immediate decay, 
0-states localize at the interfaces of large domains of spins in the states ±1. Hence the field ip(x,t) indicating the 
local presence of intermediate states is zero if 4>(x,t) is close to —1 or 1, and non-zero only when (j)(x,t) takes values 
far from the absorbing boundaries. The approximation ^(x^t) ~ 1 — <fi(x,t) 2 is thus justified and a single field <f> is 
expected to be sufficient to describe the behavior of domain walls driving the asymptotic coarsening. Within this 
approximation, and neglecting the modulation of the diffusion constant induced by the term 2 ^^^^ i n front of the 
Laplacian in the first equation of (|2|) (this term takes values near 1 close to the interfaces and 1/2 in the bulk), the 
system is described by the following equation for the local magnetization field (j>, 



—<f>(x, t) = - [4>(x, t) - <j>{x, tf] + V 2 0(x, t) + y/l-cj>(x,t) 2 r){x, t), (6) 

where rj is the usual Gaussian noise, uncorrelated both in space and time and of unit variance. The first term on the 
right-hand-side of Eq. absent in the voter model, is the typical drift term related to a double- well potential in 
the sense of a Ginzburg-Landau model. It drives the system to relax towards the two minima of the potential at ±1, 
coinciding with the absorbing states of the 3-state model. 

In the previous section we have seen that the algebraic nature of the decay of the interface density and the 
corresponding coarsening exponent are not affected significantly when the number of intermediate states is increased, 
provided that we use a model in which interfaces remain compact (e.g. when the transition rules favor the convergence 
to the absorbing states). For the multi-state model proposed, we expect Eq. ([6|) to give a valid asymptotic effective 
description in the limit of continuous spins. To this end we study the limit in which the number of intermediate states 
S is large and in which we can approximate spin configurations by a continuum profile. Specifically, consider p — jr^ji 
where the spin state s takes values s € {1, . . . , S}. In the limit S ^> 1, p can then be thought of as a continuous 
degree of freedom (at each lattice point), taking values p E [0, 1]. In order to recover a well-defined model it is then 
necessary to re-scale the spatial dimension appropriately, viz. x — *■ x/S — 1. Uniform profiles at p = or p = 1 are 
then absorbing states. In this limit the probability that the field pi at site i increases (by an infinitesimal amount) at 
any given time step is then linear in the value pj of one randomly chosen neighbor with which interaction occurs. 
The continuous representation provides a good long-time coarse-grained approximation for a model with large S 
in the coarsening regime, where L(t) is much larger than the interface width. The coarse-grained nature of this 
representation allows one now to derive a phenomcnological Langevin equation for the variable p(x,t). Considering 
a Fokker-Planck equation for a homogeneous mean-field variable p and imposing the presence of absorbing states at 
p = and p = 1 , we find 

^P(p, t) = (1 - 2p) ^ [p(l - p)P(p, t)\ + l^ W - P)P(P, t)] , (7) 

where p — J dyyP(y) is the average value and where 7 relates to the amount by which fields are incremented or 
reduced in the course of a microscopic step. The corresponding Langevin equation for the spatial dependent field 
p(x, t) then reads 



-£ = -(1 - 2p)p(l -p) + DV 2 p + Vtp(1 - p)v, (8) 



where we have added a diffusion-like term. The noise 77 is delta-correlated in time and space. The substitution 
<p(x, t) = 1 — 2p(x, t) in Eq. © yields a Langevin equation of the same form as Eq. ©. While no claim is made that 
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this derivation is rigorous or exact, these simple arguments lead one to expect that the entire class of voter models 
with S intermediate states exhibits an asymptotic coarse-grained dynamics described by Eq. ©. Since Eq. ^ leads 
to algebraic coarsening as discussed below, this is consistent with our numerical findings for the multi-state model 
(see Fig. g}. 



IV. ANALYSIS OF THE EQUATIONS 

A. Field theoretic analysis and relation to the GV class 

In this section we analyze the field theoretic counterpart of this model, that can be generally identified starting 
from the equation 

d 



— #c, t) = act>{x, t) - a<t>{x, tf + DV 2 <t>{x, t) + y/>y(l - (f>(x,t) 2 )r)(x, t), (9) 

i.e. a special case of the Langevin equation of the Generalized Voter class proposed by Al Hammal et al. in Ref. [22j 
(one has b = 0, a > in their notation). The fact that linear and cubic terms have the same coupling constant is 
necessary to enforce Z2-symmetry and has striking consequences for the corresponding field theory. 
Using the response functional technique [28[ , we obtain the following effective action 

S]d>, $]= f d d xdt 



<f>(d t -a- DV 2 ) <j) + acfxj) 3 - jcj) 2 + jcj) 2 cj) 2 . (10) 

Naive power counting reveals the canonical dimensions of the fields, [</>] = 1 and [(f)] — k d and that of coupling 
constants, [a] — k 2 and [7] = k 2 ~ d . The noise vertices become relevant for d < 2, whereas for d > 2 the process is not 
affected by multiplicative noise. In fact, the critical theory (a = 0, i.e. voter model) can be renormalized exactly to 
all orders, leading to a complete characterization of the ordering process in any dimension, above and below d c = 2 

For any a > 0, the system is in the broken symmetry phase with spontaneous magnetization and minima in ±1. In 
analogy with model A below the critical temperature, the dynamical field theory should be governed by the strong 
coupling fixed point a — > 00 [3, [25J. This can be verified by developing momentum shell renormalization for the 
coupling a. Under rescaling of lengths and times, x — * bx, t —> b z t, one considers the transformations D — ► b z ~ 2 D, 
a — > b z a and 7 — ► b z ~ d "f. Imposing scale invariance for the system's fluctuations, as required by coarsening regime, we 
get z = 2 and the tree-level flow equations ^ — 2a, and 4? = (2 — d)j, where b = e l . These simple renormalisation 
group (RG) equations hence predict a strong coupling fixed point for a, playing the same role as the zero-temperature 
fixed point does in the Ising model. This means that in the region a > we do not expect corrections to the mean- 
field dynamical exponent z = 2, so that this exponent, characterising the growth of the correlation length during 
coarsening, is expected to be accurate in any dimension. This is also true for a — since the propagator does not 
acquire corrections. The density of interfaces on the other hand follows the general scaling n(t) ~ \/L{t) ~ t~ x l z in 
any dimension only for a > 0, whereas it is modified to a logarithmic behavior for the critical (voter) theory (a = 0). 
We have here disregarded corrections to the tree level RG equations and have neglected higher-order diagrams. The 
absence of loop corrections for a > can be seen to be justified by performing a mapping of this model on branching 
annihilation random walks (BARW) by rapidity inversion [2^, [34[ ■ Our model here corresponds to unphysical B ARW 
with negative reaction rates. However, following (29j, loop corrections shift the critical value at a = to some a c < 0, 
whereas our discussion concerns the absorbing phase at a > 0. Note that region a < does not correspond to any 
microscopic model obtained adding intermediate states to the voter model, but presents a non-trivial competition 
between the two absorbing states at ±1 and the minimum of the potential at zero |22|. For this case, a complete 
understanding of the RG diagram is still lacking [2^, [H, [35| . 

The fact that noise flows to zero for d > 2 and to the strong coupling fixed point for d < 2 indicates that roughening 
phenomena may be present below the upper critical dimension d c = 2, but that they are absent in d > 2. In dimensions 
higher than two the system therefore orders like a zero-temperature Ising model. The behavior at d = 2 can be 
understood by mapping the field theory on the dynamics of an interface using standard methods [24| . Consider the 

Hamiltonian formulation for the equilibrium problem in absence of multiplicative noise Ji = J ^(V0) 2 + V(</>) d d x, 

where V{<f>) — |(1 — <j) 2 ) 2 , and define the interface position is z — h(y,t), with y being the d — 1 dimensional space 
coordinate specified by defining a preferred direction z for the interface movement (i.e. x = (y,z)). Neglecting the 
multiplicative noise, it is a zero temperature Ising model, so the interfaces will be mainly driven by diffusive forces 
[3^]. On the other hand, the multiplicative noise term behaves like additive noise at the interface (since <f> — near 
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the interface), resulting in the following Edwards- Wilkinson (EW) type of equation, 

dh(y,t) 



dt 



DV±h(y,t) + r)(y,t) 



(11) 



with (r](y,t)r](y' ,t')) = jS(t — t')5 d ~ 1 (y — y'). Therefore the interface roughening in d < 2 is described by the EW 
universality class in di n t = d — 1 dimensions. In di nt = 1, EW universality class predicts roughening exponents 
a = 1 — dint/ 2 = 1/2 and /3 = a/z — 1/4, consistent with the behavior of the interface width W measured for our 
model in simulations (see Fig. [3j. 

Interestingly, we could also apply an infinitesimal magnetic field in order to break the Z2-symmetry, then expanding 
the fields around one of the absorbing states (e.g. (f> = 4>o — <fr' with <fio = 1) 37]. In this case, the resulting action S' 
for the perturbation field 0' reads 



= / d d xdt 



<j>' (d t + 2a- L>V 2 ) <p' - 3a<j>'cj>' 2 + 2 7 0' 4>' - 7 0' ^ 



(12) 



in which we have neglected higher order terms in <j)' . 

Eq. (fl2"|) is similar to the action of directed percolation (DP) Q, therefore it is convenient to put it in the standard 
DP form [38| by rescaling the fields. Let us define £ = y/3a/2j(f)' and £ = y/2-f /3acf)' , after a bit of algebra we get 



S% C] = / d d xdt ( (d t + 2a- DV 2 ) C - u U - Q 



CC - 7C 2 C 2 



(13) 



where u — ^/Qarf. Like in DP, the pure noise coupling 7 can be neglected with respect to u, since [u] = k 2 rf / 2 and 
[7/u] = fc _d / 2 , thus the effective action becomes 



d d xdt 



((d t + 2a-DV 2 )(-u (C-C)CC 



(14) 



It is easy to check that one loop correction to the propagator shift the bare critical point (a = 0) outside the region of 
interest (i.e. towards negative values a^j < 0). In other words, when we apply a small magnetic field to our system, 
the corresponding field theory should be that of directed percolation in the absorbing or active phase depending on 
the sign of a. In terms of interface dynamics, the presence of a small magnetic field breaking the Z2-symmetry should 
bring the system in the Kardar-Parisi-Zhang universality class [24j. It would here probably be desirable to have a 
microscopic model with S > 3 states available to investigate such cases with broken Z2-symmetry further. 



B. Numerical integration 



To verify the validity of the equations proposed in the previous sections, and to compare their behavior with the 
results of microscopic simulations, we have integrated them numerically, considering both deterministic cases and 
versions of the equations which are subject to noise. It is well known that multiplicative noise prevents integration 
by naive Euler-forward discretisation, but various methods have been recently proposed to overcome this problem 
(see e.g. 35, [3!| and references therein). Unfortunately, a generalization of these step-split methods (35l . f39t ] to 
systems of coupled stochastic differential equations with cross-correlated noises is still awaited, therefore in the case 
of Eqs. ([2]) we limit our study to the deterministic variant of the equations. Slight dependencies on e.g. system size 
and time-discretisation in the numerical integration of these equations arc here hard to eliminate due to limitations 
in computing time (we have used a standard Euler-forward scheme at time stepping of 0.01). Nevertheless, the 
results for the density of interfaces n ~ 1 — (<fi 2 }, reported in Fig. [TJ are in good agreement with measurements 
from microscopic simulations, indicating a scaling exponent of about 0.45 to 0.5. They are also consistent with the 
snapshots generated from the microscopic dynamics (see Fig. [5]) in which a curvature-driven coarsening is observed, 
even if domain boundaries are not completely smooth due to the interfacial roughening discussed above. This provides 
further evidence that noise cannot be expected to play a relevant role in the asymptotic coarsening dynamics of the 
voter model with intermediate states. 

In the case of Eq. ^ the step-split method can be applied and both stochastic and deterministic versions of this 
reduced equation can hence be studied numerically. Split-step methods for the stochastic case here rely on the fact 
that the Fokker-Planck equation for stochastic processes of the form <fi = \f<\>r\ can be solved exactly, and that the 
resulting field distributions can subsequently be sampled. We here employ the scheme of [391 ] . using the replacement 
\J\ — 4> 2 — > 0(0) \/l — 4> + 0(— 0)V1 + 4> as described in [35]. The deterministic remainder of the equation is then 
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FIG. 5: Snapshots of the coarsening dynamics of the model with 5 states. Data is from simulations of a system of size 300 2 
sites, only a region of 100 x 100 spins is shown. Snapshots are taken after 200, 400, 800 and 1600 sweeps respectively (left to 
right). In the lower row we depict only spins in the intermediate states. 



integrated using a simple Euler-forward scheme. Typical time steps used are here of the order of 0.25 [35J, l39f | - 
Results are shown in Fig. [T]for different values of the noise amplitude 7. Although the quantitative details of these 
numerical findings show dependencies on parameters such as system size or time-stepping, results are consistent with 
an exponent of 0.45 to 0.5 for the decay of the density of interfaces, and hence in good agreement with the above 
theoretical considerations and direct measurements from simulations of the microscopic dynamics of the 3-state model. 
At variance with the analogous equation describing the standard voter model without intermediate states 35], and 
in agreement with the field theoretic prediction, the scaling behaviour of the coarsening process is not significantly 
affected by noise amplitude 7 in Eq. ©. In conclusion, numerical data obtained from integration of Eqs. and ^ 
corroborate the theoretical analysis and demonstrate the validity of the continuous approach to describe long time 
coarsening behavior of the microscopic voter model with intermediate states. The fact that both Eqs. ^ and Eqs. 
(I6I9|) give an exponent 6 ~ 0.45 and that this value does not change much with the amplitude of the noise, leads one 
to suggest that the discrepancy with the theoretically expected value 0.5 of Model A dynamics is probably due to 
finite size effects in the simulations, or to other numerical inaccuracies. 



V. CONCLUSIONS 

The voter model is one of the most studied microscopic models of interacting particle systems [|| , and a paradigm for 
phase-ordering processes driven by interfacial noise. Recent developments in the study of non-equilibrium statistical 
physics have put forward a common paradigm to study coarsening with and without surface tension, which has led to 
the identification of the generalized voter (GV) universality class [22J. Voter models with modified spin- flip probability 
(i.e. breaking the linearity with the local field) have been proved to belong to this general class [E US [U E3| • I n this 
paper we have shown that voter models with intermediate states fall into this class as well, and have studied their 
scaling behavior in detail. In particular, we observe algebraic domain ordering and the density of interfaces is found 
to decay with an exponent of 0.45 to 0.5. 

We have derived a set of two coupled stochastic differential equations for appropriately defined fields, and have demon- 
strated that these equations faithfully describe the long-time large-scale behavior of the 3-state voter model. Within 
a physically motivated approximation the two equations can be seen to reduce to one equation for a single field, and 
the latter equation in turn takes a form of a special case of the GV class. Moreover, the same algebraic coarsening 
is found in models with more than one intermediate state, and similar equations appear to hold for such cases with 
S > 3. 

We have used a field theoretical approach to understand the universal properties of the model and have demonstrated 
that numerical results (obtained from both microscopic simulations and numerical integration of the effective equa- 
tions) are consistent with this theoretical analysis. Interestingly we find that, due to the appearance of an effective 
surface tension induced by the presence of intermediate states, the domains coarsen like in the low temperature phase 
of Model A dynamics. In d < 2, the presence of interfacial noise induce interfacial roughening in agreement with 
usual Edwards- Wilkinson scaling. RG equations predict smooth interfaces in higher dimension. 

It is hoped that our theoretical analysis will help to shed light on the phenomenology of this class of generalised 
voter models. These models have come to the attention of statistical physicists only recently, mostly due to their 
applications in the modelling of social systems. With a modest amount of imagination however, they can be expected 
to be applicable also in the study of kinetic reactions with intermediate unstable compounds and in other physical, 
biological or chemical non-equilibrium systems. 
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